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Abstract 

We formally state the connection between the relativistic part of the radiation reaction and the Poynt- 
ing Robertson force term, —Rv/c 2 , where R is power radiated. Then we address the question, does 
[x,p] = ih for an accelerating charge ? The full radiation reaction term is used, which includes the 
relativistic term (von Laue vector.). We show that the full relativistic radiation reaction term must be 
taken into account if a commutation relation between x and p is to hold for an electron under uniform 
acceleration, consistent with the expectation values of x 2 and p 2 . 

PACS codes: 41.60.-m, 41.75.Ht, 41.20.-q, 03.65.-w 

Key Words: [ Commutation relations, accelerating electron, radiation reaction, von Laue T, relativistic 
electron, LAD equation, Poynting Robertson drag.] 



Introduction 

It has been shown [1] , that the nonrelativistic theory of the electron is fundamentally inconsistent unless 
both radiation reaction and the vacuum field are allowed to act on the electron. In particular, it was 
shown that the canonical commutation rule [x,p] = ih for the electron is violated if radiation reaction 
is ignored. In this previous work [1], the equation of motion used was that of Lorentz. 

r(t)-ir(t) = -E (t) (1) 
m 

where 7 = 2e 2 /3mc 3 and r is a vector. We instead use the relativistic LAD equation of motion (of 
Lorentz, Abraham and Dirac) [2, 3, 4] given by 

x v - j(x v + x^Xftiv/c 2 ) = —x^Fp (2) 

where the damping term 7 = 2ro/3c is defined above and has dimensions of time, and rg = e 2 /mc is 
the classical electron radius. The term on the right is equivalent to eE(t)/m where E(t) is the electric 



field. Also m is the rest mass, which we shall use throughout the paper. 

Let us examine this equation (2) a little more closely. The first term clearly comes from the kinetic 
energy of the charged particle (electron in this case) . The second term on the left is the Abraham Lorentz 
radiation reaction term F = 2e 2 d/3c 3 or it is derived from what used to be called the acceleration energy 
by Schott 1915 [5]. The third term on the left is of the form F = -Rv/c 2 where R = 2e 2 a 2 /3c 3 which 
is the regular Larmor power radiated by an electron of acceleration a. This term corresponds exactly 
to the Poynting-Robertson force term [6, 7]. This is exactly the Abraham result of 1903 [3] and the 
relativistic part of the von Laue 1909 radiation reaction T [8]. A history of the Dirac equation (2) is 
given by Rohrlich [9] and the book by Grandy [10]. The remarkable thing about Grandy's book is that 
on page 204 he writes the drag term (in the LAD equation) in the format of the Poynting-Robertson 
drag —Rv/c 2 , where R is the power radiated by the accelerating charge. However, no connection is 
made between the (relativistic) radiation reaction term and the Poynting-Robertson drag force. The 
author feels that this is an oversight, no electromagnetic text book mentions the connection even though 
the relativistic term is very important and used all the time, possibly without people realising it. The 
two radiation reaction force terms are both derived from the same power expression. Usually text 
books (for example Jackson [11] or Griffiths [12]) do not derive force expressions. They typically start 
with the electric field of a moving charge, then calculate the Poynting vector and then integrate over a 
cross-sectional area to find the power radiated into a given angle. Very few text books even mention the 
famous Dirac paper of 1938. Jackson's book refers to it in the problems section only, Griffiths' book does 
not mention it at all. It should also be mentioned that Dirac's 1938 paper contains the original discus- 
sion of the advanced and retarded fields, which was elaborated on by Wheeler and Feynman later in 1945. 

When acceleration is taken into account, it is most important to keep the relativistic radiation reaction 
term. Boyer [13], noted that an electric dipole accelerated through the vacuum would see a surrounding 
field not quite equal to the usual Planck distribution. A correction term was needed [14], which turned 
out to be exactly the relativistic radiation reaction term, (or the Poynting and Robertson drag). Boyer 
[14], showed that a 

" classical electric dipole oscillator accelerating though classical electromagnetic zero-point 
radiation responds just as would a dipole oscillator in an inertial frame in a classical thermal 
radiation with Planck's spectrum at temperature T = ha/2Trck" 

where T is the Unruh-Davies temperature. 
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Later, Boyer [15], did a similar calculation for the spinning magnetic dipole and found a mismatch with 
the Planck distribution again. He later corrected the magnetic dipole work with a similar drag force to 
regain the Planck distribution, [16]. He was able to show that, 

" the departure from Planckian form is canceled by additional terms arising in the relativis- 
tic radiative damping for the accelerating dipole. Thus the accelerating dipole behaves at 
equilibrium as though in an inertial frame bathed by exactly Planck's spectrum including 
zero-point radiation." 

We have established that for an accelerated charge, consistency with thermodynamics requires that we 
keep the relativistic radiation reaction term in the equation of motion of the electron. We now turn 
our attention to the [x,p] calculation. Instead of the quantized vacuum field, we will be dealing with 
the vacuum field and the Planck distribution at a temperature given by the Unruh Davies temperature 
T = ha/(2nck B ) where a is the uniform acceleration and k B is Boltzmann's constant [14, 17]. 

We are strongly motivated by Wigner [18], who once asked the question; 

"Do the equations of motion determine the quantum mechanical commutation relations? " 

Wigner showed that the equations of motion for a harmonic oscillator do not uniquely determine the 
commutation relations. We find that indeed the equation of motion in our case does determine the 
commutation relation for x and p. The commutation relation found is consistent with the expectation 
values of x 2 and p 2 but at first sight it appears it is not [x,p\ = ifi. Our commutation relation does 
however reduce to this expected result for the case when the acceleration of the charge a — > 0. 

Results and Discussion 

Consider an electron undergoing an oscillation, like a mass on a spring, with the motion restricted to the 
x direction. We shall take the resonant frequency of this oscillation to be cj - This oscillatory motion is 
taken to be nonrelativistic and hence we continue to use the rest mass of the electron m. We consider 
the amplitude of this motion to be very small (as in the small dipole approximation) and we only keep 
terms linear in the amplitude x of the oscillation. The electron is accelerating with uniform acceleration 
a in the z direction. We may rewrite the equation of motion in this case as, 

2 

x + wlx-fix - %x) = — E(t) (3) 
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where E(t) is taken to be the electric field at the equilibrium point of the oscillation, instantaneously 
at rest with respect to the oscillation, as in Ref. [14]. Note that we have ignored terms in the radiative 
damping that go like 'x 2 x since these would involve the amplitude of the oscillation squared which is 
a small quantity. We have only kept the z-component of the acceleration squared term because this 
term is only linear in the oscillation amplitude. We will write the electric field, (in the rest frame of the 
equilibrium point of the oscillation) , in quantized form as 



E(t) 



V* Jo 



du e ks (e(uj k )a ks e tUJkt + e* (w k )rf ks e lUJkt ) 



(4) 



where a ks is the photon annihilation operator for mode (k, s) and e ks is the polarization unit vector. 
The Fourier components e(uj) are well defined by Boyer [14] and in earlier work by the same author cited 
therein. We leave the definition until later. Note also that we have taken and not 1/\/2tt since 

that would imply an integral from — oo to oo. Here we wish to keep the positive and negative frequency 
components separate. Let 

x(t) = Xl e- lUJkt + x 2 e lUJkt (5) 

and treat the positive and negative frequency terms separately in the equation of motion. Substituting 
xi into our Eq.(3) gives 

e ks e(uj k )a ks 



similarly for xi we find 



x\ = - 



x 2 = - 



duj k 



(uj 2 - wg) + ij(oj s k + a 2 uj k /c 2 ) 

e fes£*( w fe) a1 ks 



thus 



x(t) 



dui k e ks 



o 



K - w o) - h(u k + a 2 uj k /c 2 ) 
e(uj k )a ks e- i0Jkt e*(uj k )ah s e iuJkt 



(^1 - w o) + + a 2 u k /c 2 ) (uj 2 - uj 2 ) - i 7 (w| + a 2 uj k /c 2 ) 



(6) 



(7) 



(8) 



From this result we can easily find 

e 



x(t) = i 



sprim J 



duj' k e k i s iuj' k 



e(iv' k )a k ' s 'e 



-iU>U 



- y^k) - k's'Z k 



K 2 - wo) + *7K 3 + ^uj'Jc 2 ) (uj 12 - uj 2 ) - i^ujf + a 2 uj' k /c 2 ) 



(9) 



The commutator then becomes 



[x,p] 



[x, mx\ 
ie 2 



duii 



duj' k e ks ■ e kls ,£(u] k )£*(ui' k )uj' k 
[a ks ,ah, s ,]e-^ k -^ 



[(uj 2 uj 2 ) + i 7 ( w 3 + a 2 uj k /c 2 )][(uj' 2 uj 2 ) i 7 ( W * + a 2 uj' k /c 2 )] 
[ot fca ,o fc / a /]e i ( w *- w i) t 



[{uj 2 uj 2 ) i^l + a 2 uj k /c 2 )}[(uj' 2 uj 2 ) + ?7 K 3 + a 2 JJc 2 



(10) 
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We have assumed here it is sufficient to use, [eifc s , ' y v] = ^ss'Sf, k , and [a^ks, a>k's'] = ~^ss'Sf.y we find 
by setting uj k = u> (we drop the k subscript which is no longer needed), 

2uj{s{uj)s*{uj)) 



ie 2 f°° 

[x,p] = / dto 

irm J 



(w 2 -wg) 2 +7 2 w 6 [l + (a/cw) 2 ] 2 
At this point we introduce the expectation of the fields given by Boyer [14], in Eq.(25) of that paper. 

ATiuj 3 



(11) 



(e(u>)e*(u>)) 



3c 3 
8^ 
3c 





coth 


' 7TC0J ' 




\OjjJ 




. a . 



(-) 

\Cld/ 



p(oj) coth 



■KCIjJ 



a 



(12) 



where p(cj) — Huj 3 / (2tt 2 c 2 ) the energy density of free space. It is important to note at this stage that 



- coth 
2 







(e 2vcUJ ' a + 1\ 


. a . 


\ 


i e 27TcL0/a — I J 



(13) 



2 e 27TcLU/a _ I 

where 27rcw/a = tku/(k B T) which implies that T — ha/ (^Trck-s) which is the Unruh Davies tempera- 
ture. The p(cj) coth term is therefore the vacuum field plus the Planck distribution at temperature T. 
Substituting the Eq.(12) into the integrand of Eq.(ll) the commutator becomes, 



M = ^ j 

K Jo 



duj 



i + coth[2fi 



(14) 



(c 2 -oj 2 ) 2 + ^ [i + (&y 

Then using the substitution z = uj 2 — cj 2 , and s = a/ (ujqc) and setting lo — oj everywhere else we get, 



[x,p] 



where A = 7 Wq(1 + s 2 ). Hence, 



2ih f°° , 7^(1 + s 2 ) coth [tt/s] 
dz- 



2ih r 

^ Jo 

2ih , . . . r 

= coth [n/s\ I 

n Jo 



z 2 + 7 2 w6(l + s 2 ) 2 
A 



dz 



z 2 + A 2 



(15) 
(16) 



[x,p] = — coth [tt/s] [tan 1 (z/A)]f° 



7T 
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ih coth [TTcaJn/a] 



(17) 



where we have substituted back for s in the last step. This is not quite the expected result, but we do 
know that when a — > then coth[oo] = 1 and we get back the expected result. To show that this is 
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indeed the correct result we need to consider the uncertainty relation, AxAp. Using, 

Ax = yV) - (x) 2 

Ap = m\J (x 2 ) — (x) 2 

(18) 

and (x) = (x) = 0, using the results from Boyer [14] we find; 



AxAp > m\J (x 2 )(x 2 ) 

> ^ coth [ncwo/a] (19) 

Generally speaking, we know from quantum mechanics text books, that if we have the uncertainty rela- 
tion AaAb > \d\/2 then this implies a commutator [a, b] = id. In our case above d — ?icoth [ncojo/a]. 



Conclusions 

In conclusion, we have made the connection between the relativistic radiation reaction term (in the LAD 
equation of motion) and the Poynting-Robertson drag force term. We have further shown that it is 
necessary to use the full equation of motion given by Lorentz, Abraham and Dirac (LAD) when treating 
an accelerated electron in order that the commutation relation (which agrees with the corresponding 
uncertainty relation) between x and p hold. The result reduces to the expected result when the acceler- 
ation a — > 0. 

We also note that we have assumed a form of electric field quantization without a rigorous deriva- 
tion. Our end result does appear consistent with previous work. However, we note that to regain the 
regular commutation relation it may be necessary to consider the relativistic Dirac equation, and use a 
combination of a or 7 matrices rather than x and p operators. 

Finally, we note that is it possible to derive the result quickly by using the vacuum thermofield oper- 
ators originally defined by Takahashi and Umezawa [19, 20] for the a and a) operators which then gives 
the coth [-kcwq/cl] term directly in the commutator. The thermofield commutator for [a,k s , oVs']t = 
coth [ttcloq/ a] 5 ss >5l k , would replace the commutator above Eq.(ll), where we have used the Unruh 
Davies temperature T. 

The Thermofield operators use a vacuum double Hilbert space, with a fictitious mode or virtual photon 
mode for the vacuum. The ordinary a and a* operators are transformed via a Bogoliubov transformation 
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from a unitary operator T(9). The thermofield vacuum is given by, 

|0> T = T(0)|OO> (20) 

where T{9) = cxp[— 9(aa — a^a* 1 )}, where the tilde denotes a fictitious or virtual mode. Note that 
Tt(0) = T(-6). Using 

e A Be- A = B + [A,B] + ±[A, [A, B]] + i[4 [A, [A, B]]] + (21) 

we find, 

a T = T{9)aT jf (9) = a cosh 9 - ^ sinh 9 
a} T = T(0)atT\6) = a f cosh 9 - a sinh 9. (22) 

(23) 

It is easy to show that the photon number and commutator for the thermofield operators become, 

(00V T a T |00) = sinh 2 9 = — ^ - 
(00|[a T , a) T ] |00) = sinh 2 9 + cosh 2 9 = coth a (24) 

(25) 

where a = TTUj c/a. Clearly we could have used these for our accelerating frame operators and saved 
a bit of time calculating electric fields like Boycr [14]. It is not clear if the thermofield operation can 
be related to space-time curvature, and quantum gravity, but the accelerated frame appears to give the 
same result as would the thermal frame using thermofield operators and the Unruh Davies temperature 
relation. It may be that we should simply normalize the [a, a'] commutator to account for the thermal 
result, coth [ncujo/a], in which case you regain the usual [x,p] commutator. This would also normalize 
the uncertainty relation and we would get back AxAp > Ti/2. 
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